Graphene with Structure-Induced Spin-Orbit Coupling: 
Spin-Polarized States, Spin Zero Modes, and Quantum Hall Effect 
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Spin splitting of the energy spectrum of single-layer graphene on Au/Ni(lll) substrate has been 
recently reported. I show that eigenstates of spin-orbit coupled graphene are polarized in-plane 
and perpendicular to electron momentum k; the magnitude of spin polarization S vanishes when 
k — > 0. In a perpendicular magnetic field B, S is parallel to B, and two zero modes emerge in the 
Landau level spectrum. Singular B-dependence of their magnetization suggests existence of a novel 
magnetic instability. They also manifest themselves in a new unconventional quantum Hall effect. 

PACS numbers: 71.7D.Ej, 73.43.-f, 73.61. Wp, 75.70.Cn 



Since the discovery of graphene with its quasi- 
relativistic energy spectrum of zero-gap Dirac fermions 
[3, 0] and unconventional Quantum Hall Effect (QHE) 
in single- and bilayer graphene 0, 0, @, this mate- 
rial attracts attention because of its unique electronic 
properties and prospective applications in nanoelectron- 
ics p], 0|. Applications in spintronics depend strongly 
on the control of spin-orbit coupling. Spin-orbit cou- 
pling in graphene comprises intrinsic and extrinsic com- 
ponents [1, 0]. First one is known to be very weak 
in plane graphene [Tot ITU; in carbon nanotubes it is 
due to their curvature [12 . Il3 |. Extrinsic spin-orbit 
coupling originates from the interface between graphene 



14] 



and substrate. Recently Varykhalov et al. 
ported spin angle-resolved photoemission spectroscopy 
(SARPES) data taken from single-layer graphene on 
Ni(lll) intercalated with a Au monolayer. They re- 
veal strong momentum-dependent in-plane spin polar- 
ization. The monolayer of Au maintains the symme- 
try of graphene and strongly reduces the threefold de- 
formation of graphene (originating from its coupling to 
Ni substrate) and nonequivalence of MB) lattice sites 
resulting from it. According to Ref. llj, the technique 
brings the system closer to ideal freestanding graphene 
than any other preparation on a solid substrate before. 
These results call for a systematic theory of spin polar- 
ization in freestanding graphene with extrinsic spin-orbit 
coupling, and in this paper I provide such a theory for 
free electrons and electrons in a quantizing magnetic field 
B. The results are in general agreement with the data 



by Varykhalov et al. [14| and predict (i) dependence of 



spin polarization on the magnitude of the momentum k 
at the spin-orbit momentum scale, (ii) spin zero modes in 
a quantizing field indicating a novel magnetic instability, 
and (iii) a new unconventional QHE. 



A 4 x 4 Hamiltonian of graphene with extrinsic spin- 
orbit coupling can be conveniently represented in terms 



of a Kronecker product of 2 x 2 matrices cr and s as Q 



H° K = T (<7 • k) + ~A(<7 X s) 



(1) 



Here 7 = hi 



10 8 cm s 1 , and A is a spin-orbit cou- 



pling constant; for graphcne/Au/Ni(lll), A « 13 meV 
jl4| . Pauli matrices of pseudospin cr — (<xi,<72) operate 
on A(B) lattice cites, and s = (s x , s y ) are Pauli matrices 
for real spin. It is seen from first term in |T]) that cr trans- 
forms as a momentum k. Therefore, as was mentioned 
by Kane and Mele Q, spin-orbit term can be considered 
as a modified Rashba term with k — > cr; the conventional 
Rashba term (k x s) is small and will be disregarded. In 
the same representation, 4x4 spin matrices are cos, with 
<7o for a unit matrix in cr space. The subscript K in ri° K 
indicates that it is acting in the vicinity of K point of the 
Brillouin zone, FIG. 1; for K' point, H° K , = aiV^a^ 1 . 

It is convenient to change from the 4x4toa2x2 
formalism. First, we perform an unitary transformation 
of the Hamiltonian, TCk — Uri^U" 1 , and of spin matri- 
ces <7o s with an unitary matrix U = [(sq + s z )ctq + (sq — 
Sz)< J i]/2, sq being a unit matrix in spin space. Then 
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c 

As, 



S = (Ti s, S z 



(2) 



with C = "f(sok x — is z k y ); here (S,^) are new spin 
matrices. Next step is eliminating lower components of 



eigen-bispinors ^ = 



(e-Asj,) -1 ^^ 



of the equation H^f = £\1/ 



7 



ek 4 



-iAfc_ 



e 2 — A 2 V i\k + efc_ 



Then we arrive at a 2 x 2 Hamiltonian 

, . _ 7 2 / ek-k + -i\k 2 _ 
H ^ £ > ~ £ 2„ A 2 I i\k 2 ek+k^ 



(3) 



(4) 



that depends explicitly on its eigenvalues e, H(e)<p = eip. 
When B ^ 0, operators k± — k x ± ik y do not commute; 



2 



below, B = Bz. We note that Eq. ([4]) is exact, and 
despite the fact that the problem is nonlinear in e, calcu- 
lation of spin polarization in this representation is more 
concise than in the original 4x4 representation. 
For free carriers, B = 0, eigenspinors of H(e) are 



1 

V2 



ivk 2 _jk 2 



u = ±l. 



(5) 



From © and Q follows an equation e v (e v + \v) = -f 2 k 2 
for the eigenvalues e v . Its solutions are 



M*0 = V - (-A + M VA 2 + 4 7 2 A; 2 ) , /i 



±1. 



(6) 



The spectrum includes two zero-gap branches and two 
gapped branches of the same shape. The gap equals 
2 1 A | , and the separation between gapped and ungapped 
branches is fc-independent and equals A. The spectrum is 
the same as for unbiased bilayer graphene without spin- 
orbit coupling [IB], [lfl Il7l |. but the region of parabolic 
expansion, k <C |A|/7, is narrow because A is small. 

It is easily seen from ([2]) and ([5|) that {<p v \S z \ip u ) = 0, 
hence, spins are in-plane polarized. Because S includes 
o i , calculation of in-plane polarization involves the lower 
spinor (f> and is more cumbersome. Nevertheless, it is 
straightforward, and applying (J3]), ([5]) and ([6]) results in 
in-plane spin polarizations S 1 , Al (k) for all (v, jj) states 



(k) 



<*^|S|*,M> 



2/x7(k x z) 
(IVl^V) ~ \A 2 + 4 7 2fc2 



(7) 



Eq. ([7]) indicates transverse spin polarization, FIG. 1, as 
concluded by Varykhalov et al. 14]. Its magnitude is 
fc-dependent. When k ^S> k\, k\ = |A|/2 7 being a char- 
acteristic spin-orbit momentum, it saturates, |S„ M | — ► 1. 
In the opposite limit, k <C k\, it vanishes as k/k\. 

Chirality of the spinor </?„(k) is defined by v, spin po- 
larization S„ M by /I, and the product fii/ specifies elec- 
tron and hole spectrum branches. Experiments of Ref.ll4 
measured the magnitude of A, |A| ~ 13 meV. Measuring 
the sign of would allow finding the sign of A; indeed, 
it is seen from (J6j) that /i/A > for external Fermi circles. 
Due to the requirements of time-inversion symmetry, spin 
polarization is identical near K and K' points. It is not 
clear currently which of the factors (experimental resolu- 
tion, temperature, or staggered potential of Ni substrate) 
was the main obstacle for measuring SARPES spectra for 
k < k\. However, measurment of S ufl for k > k\ should 
shed additional light on the role of these factors. 

Application of well developed techniques for detecting 
in-plane polarization S(k), based on Kerr spectroscopy 
[l8L Hil . [2~oj ] and spin-galvanic effect 21], is hampered 
by the conductivity of metallic substrate. Reducing the 
thickness of the substrate to only a few monolayers or 
developing insulating substrates can render them proper 
efficacy. 



For B || z, applying a Peierls substitution k = — iV + 
eA./hc, A being a vector potential, one expresses k± in 
terms of Bose operators, k + = (^/2/£)a + , fc_ = (^/2/£)a, 
[a,a + ] = 1; here £ = -J 1 chjeB is a magnetic length. 
Then, instead of ([3]), one arrives at 



H(e) = |A| 



2T 2 



-i(3a 2 



i(3(a+) 2 



(8) 



where (3 = A/|A|, e = e/|A|, and T = j/£\X\. The Hamil- 
tonian H (e) depends explicitly on its eigenvalues e. From 
here on, energy e is measured in the units of |A|. 

Solution of the corresponding eigenspinor problem can 
be found in terms of oscillator eigenfunctions tp n p3 | 



Cllpn- 

rV2 f t\jn — 



1 



1 ci - ifi^fn c 2 
- 1 ci + ei/n c 2 



(9) 



for n > 2. The coefficients ci,2 = ci^(n), normalized as 
|ci| 2 + |c 2 | 2 = 1, read as 



ci i/3e(l + 2nT 2 - e 2 ) 



-, c 2 = 



(10) 



c 2 2y/n(n-l)T 2 ' 2n(l + r 2 )-e 2 

Eigenvalues obey the equation 

e 4 - [1 + 2r 2 (2n - l)]e 2 + 4n(n - l)r 4 = (11 



and are 



(O 2 - 



1 



1 



4r 4 



2(2n - l)r 2 ± a/1 + 4(2n- 1)T 2 

(12) 

Eq. (|12p coincides with the expression for bilayer 
graphene in absence of spin-orbit coupling 17 , 23|, [2^ . 

In addition to the solutions with n > 2, there are so- 
lutions with n < 1 that are of special interest because 
of their peculiar spin properties. For n = 0, there is a 
single solution ci(0) = 0, c 2 (0) = 1, (f>o = 0, with e = 0. 
For n = 1, there are three solutions with ci(l) = and 
c 2 (l) = 1, they differ by their (f>i spinors. For one of them 
eigenvalue vanishes, e° = 0, and components of the spinor 
4>1 are C3 = ipr^/2, 04 = 0. Two nonvanishing eigenval- 
ues are ei = ±Vl + 2r 2 ; the components of the corre- 
sponding spinors <f>i are C3 = —i/3/T\/2, C4 = ei/T^/2. 
These expressions can be also found from ([9]) and from 
(fT2| [with the upper sign in (fT2l ] by plugging n = 1. 
Hence, there exist two zero modes, eo = and = 0, 
and spin-orbit coupling of (er x s)-type preserves two-fold 
degeneracy of e — state typical of single-layer graphene 
[EiljIUlj but changes its nature: Degeneracy is dynam- 
ical rather than Kramers (small Zeeman splitting [26[ is 
disregarded) . 

It follows from (|9|) that in-plane spin polarization S 
vanishes in all eigenstates. Indeed, because of the fac- 
tor <j\ in S = ctiS, it mixes different components of the 
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bispinor "J", ip and <fr, and mean value of their product is 
proportional to scalar products of the oscillator functions 
tj) m with quantum numbers that never coincide. There- 
fore, S„ = for all n. This result is expected because 
Sj/ M (k) of (J7J) vanishes after averaging over the direction 
of k, and follows from axial symmetry of the problem. 

Longitudinal polarization S z , after eliminating (f> com- 
ponent of can be expressed in terms of its <p component 

S z = (<f s z + -^—j[s Q + (aa + + a + a)s z } <p\ 

(13) 

where e is energy of the eigenstate cp. The normalization 
factor equals 



(#|#) = 1 



20[( 



(e 2 - l) 2 
+ (e 2 + 1) [(aa + + a + a)s + 



2 2 1 

a ' ) s_ — a s+ 1 



(14) 



Explicit expressions for (| 13j) and (JT4J) follow from (J9J) and 
()10jl . Their original form is cumbersome but greatly sim- 
plifies after higher powers of e 2 are eliminated by employ- 
ing (jTTJ) , and (3 2 = 1 is applied. The final equation, when 
expressed in terms of e^, is rather concise 



2r 2 (2nr 2 - e 2 ) 



e 2 (l-2r 2 ) + 2nr 2 (l + 2r 2 ) 



(15) 



This equation, together with (|12|) . provides exact expres- 
sions for independence of spin polarization for all states 
with n > 1 and e ^= 0. Because e appears in (|15p only 
as e 2 , polarization is charge symmetrical: It coincides 
for electron and hole states with the same n and le^l. 



In the weak field limit, 
t ~ 1 and e, 



2^7 



field limit, T 2 



ttT 2 <C 1 , e 2 , eigenvalues are 
l)r 2 , and S z is n- independent 
and proportional to B, (Sf)^ w =F2r 2 (n > 1 for 
e~ states). Because of the spectrum degeneracy at 
k = point, the sequence e~ is nonequidistant despite 
the parabolicity of the spectrum [22}. In the strong 
► 00, eigenvalues are e~ ~ Yy/2n and 
n — 1), and spin magnetization saturates, 
(S z )n ~ =Fl. In this limit, two ladders nearly overlap and 
are split by Ae„ s» y / n72/(2r); the splitting increases 
with n but for T 2 3> n is small compared with the level 
separation Y j \pl~n inside each ladder. For A « 13 meV, 
the field separating these two limit cases, found from the 
condition Y = 1, equals B CT = c\ 2 /{ehv 2 ) » 0.3 T. 

In both limits, spin magnetization has opposite sign for 
e + and e~ ladders, hence, magnetization oscillates when 
Landau levels cross the Fermi level. These de Haas-van 
Alphen type oscillations can be detected by Kerr effect 
spectroscopy |27| even for ferromagnetic substrates. 

Semiclassical regime is achieved for Y <C 1 and n ^> 1 
with nr 2 = k 2 /2 = const; for k = 7fc/|A|, one recoveres 
© from (H2J). Keeping k = const, one finds for electron 



branches, > 0, in the first order in 1/n, 
et « - (Vi + 4k 2 ±1 



{S z )t « T 



2nVl + 4k 2 ' 

k 2 _ 2r 2 

iVl + 4k 2 ~ T a/1 + 4( 7 A:/A) : 



(16) 



Therefore, in the leading order of the expansion, level 
splitting Ae„ « 1 remains n- independent, while spin po- 
larization S z has opposite sign for two spectrum branches 
and decreases with n (or, for Y= const, with the electron 
momentum k). The oscillatory dependence of the total 
spin magnetization on B can be detected through Kerr 
spectra as discussed in the previous paragraph. 

Therefore, all quantum states with ^ are nonde- 
generate, and only two states eo = e i = are degenerate. 
These two states differ strongly in their spin magneti- 
zation. For the n = state, polarization (S z )q = — 1. 
It does not depend on B, but changes abruptly when B 
changes sign. For the n = 1 state with e\ = 0, polar- 
ization equals (S z )i = -(1 - 2r 2 )/(l + 2r 2 ); it tends to 
— 1 for B — > 0, changes sign at Y 2 = 1/2, and saturates 
to +1 for B — * 00. In this limit, contributions of two 
e = states cancel. Such a behavior is unique because 
it suggests that spin-orbit coupled graphene with filled 
e = states is unstable to ferromagnetic ordering in z 
direction in weak fields B || z, while with increasing B 
the magnetization gradually vanishes. This magnetiza- 
tion differs drastically from the edge-state magnetization 
of graphene zigzag nanoribbons proposed by Fujita et al. 
[281 ] because it originates from spin-orbit coupling rather 
than from exchange interaction. It also bears no similar- 
ity with the Dzyaloshinskii-Moriya weak ferromagnetism 
[23 . [30j ] because it is dynamical rather than symmetry 
conditioned, and develops in a paramagnet without any 
magnetic structure. Effect of electron-electron interac- 
tion on this peculiar state needs a special investigation. 

For comparison with experimental data, one needs to 
add magnetization S z of the electrons in K and K' val- 
leys. Because H° K , is related to HP K by a o\ canonical 
transformation leaving the operator S z of |[2]) unchanged, 
magnetization has the same magnitude and sign in both 
valleys. This can be also inferred from the fact that K 
and K' valleys are related by a 7r/3 rotation about the 
z axis that does not change S z . Therefore, the total S z 
equals the magnetization of K valley multiplied by a fac- 
tor of 2. 

Novoselov et al. [5j compared the conventional QHE 
[IH and two types of unconventional QHEs typical of 
single-layer 0, [j] and bilayer graphene, the material no- 



torious for its exotic QHE properties [32|. Spin-orbit 
coupled single-layer graphene introduces one more type 
of unconventional QHE. Due to the four- fold degener- 
acy of zero modes [eo = = degeneracy times fac- 
tor of 2 from isospin, K(K') valleys] the step in a xy at 
B = equals 4e 2 //i, as in single-layer graphene without 
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spin-orbit coupling. However, because spin degeneracy 
is lifted in each valley by spin-orbit interaction, and only 
isospin degeneracy persists, all B ^ steps are of 2e 2 /h. 
Therefore, the ratio of the magnitudes of B = and 
B 7^ steps equals 2, as in bilayer graphene without 
spin-orbit coupling. Depending on the magnitude of A, 
these 2e 2 /h steps can appear in pairs, as resolved 4e 2 //i 
steps, similarly to the resolution of two spin components 
of the traditional QHE. Spin-orbit coupled graphene on 
an isolating substrate would become the optimal object 
for observing this new QHE and for spin manipulation, 
a challenging task for semiconductor spintronics. 

Apparently, among the perturbations that lower the 
symmetry of the Hamiltonian the largest is the staggering 
sublattice potential of the substrate violating the equiv- 
alence of A(B) lattice cites; manifestation of the broken 
six-fold symmetry was reported in Ref. IbJ. A staggered 
potential can be described by a term Ti s t = u<j z t z , where 
t z is a Pauli matrix in the isospin space ||. This term 
creates a gap in the spectrum and lifts the symmetry of 
K(K') valleys. Adding H B t to H a K ^ K ,-, does not change 
wave functions of n — modes but changes their ener- 
gies to e — ±u (in dimensional units). For the n = 1 
soft mode, eigenvectors also change; now c\ = but 
C2,C3,C4 7^ 0. Energy spectrum can be found from a 
cubic equation that for |u/A| <c 1 defines soft modes 
£i « ±w(l - 2r 2 )/(l + 2r 2 ). As a result, the Ae 2 /h step 
in a xy splits into a plateau near B = and two e 2 /h 
steps on both sides of it. Lifting the K(K') degeneracy 
also splits all 2e 2 /h steps, making the QHE of spin-orbit 
coupled graphene similar to the traditional QHE with a 
resolved Zeeman splitting. We notice that A = 13 meV 
corresponds to a field B = 120 T for a Lande factor g = 2. 

In conclusion, a theory of the energy spectrum and 
spin polarization in single-layer graphene a subject to 
a substrate-induced spin-orbit coupling is presented. 
Energy spectrum consists of two zero-gap bands and 
two gapped bands reminding the spectrum of bilayer 
graphene without spin-orbit coupling. However, all states 
are in-plane spin-polarized perpendicular to the momen- 
tum k. This polarization saturates at large k and van- 
ishes at the scale of spin-orbit energy when k — ► 0. In a 
perpendicular magnetic field, two zero modes develop in 
each of K(K') valleys. These modes show a peculiar mag- 
netic behavior suggesting a possibility of a perpendicular- 
to-plane spin-orbit conditioned magnetism, and produce 
a new unconventional quantum Hall effect. 

I am grateful to F. Kuemmeth, C. M. Marcus, and J. 
R. Williams for stimulating discussions. 
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FIG. 1. Spin polarization of energy spectrum of spin- 
orbit coupled single-layer graphene with the dispersion 
law of Eq. © and spin-orbit coupling constant A > 0; 
energy e > A. Left column - K' valley, right column - 
K valley. Upper row - electrons, [iv = +1; lower row 
- holes, [iv = —1. External circles - fi — +1; internal 
circles - [i — —1. For A < 0, quantum numbers, [i and v, 
and spin polarizations change their signs for all branches. 
Brillouin zone of graphene is also shown. 



